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Abstract 

Functional networks, i.e. networks representing dynamic relationships be¬ 
tween the components of a complex system, have been instrumental for our 
understanding of, among others, the human brain. Due to limited data avail¬ 
ability, the multi-layer nature of numerous functional networks has hitherto 
been neglected, and nodes are endowed with a single type of links even when 
multiple relationships coexist at different physical levels. A relevant problem 
is the assessment of the benehts yielded by studying a multi-layer functional 
network, against the simplicity guaranteed by the reconstruction and use of 
the corresponding single layer projection. Here, I tackle this issue by using as 
a test case, the functional network representing the dynamics of delay prop¬ 
agation through European airports. Neglecting the multi-layer structure of 
a functional network has dramatic consequences on our understanding of the 
underlying system, a fact to be taken into account when a projection is the 
only available information. 

Keywords: Complex networks, functional networks, multi-layer networks, 
air transport 


1. Introduction 

In the early stages of complex network theory DEI. such paradigm was 
mainly used to analyze systems whose structure, either physical or virtual, 
could be directly mapped into a network: examples include transportation 
systems Pi, social networks |5] or the World Wide Web [6]. It was soon 
clear that in certain cases this was not possible, as the only information 
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obtainable from the system itself was the evolution through time of some 
observables. Such measurable variables reflect the behavior of the interacting 
elements constituting the system, and as such, the value of every observable 
is expected to be a function of the values of other peers. When the structure 
of such interactions is inferred from the dynamics of the observables, the 
result is then called a functional network. The introduction of this latter 
representation has resulted in an important step forward in network science, 
allowing a broader focus including both structural and dynamical (functional) 
relations. 

Among the examples of functional network representations, the study of 
the brain mi and of gene expressions PHU] are probably the most paradig¬ 
matic. In the former, nodes correspond to sensors of a machine recording the 
activity of the brain. Through the magnetic or electric held generated by 
spiking neurons, links are then established whenever some kind of synchro¬ 
nization is detected between the recorded time series, usually by means of 
metrics like Pearson’s linear correlation. Synchronization Likelihood HU , or 
Granger Causality na. Similarly, nodes in gene co-expression networks rep¬ 
resent individual genes, pairwise connected when some correlation is detected 
in the dynamics of their respective expression levels. 

In the last few years researchers have realized that interactions between 
the constituting elements of complex systems seldom develop on a single 
channel. For instance, in a social network, information exchange may happen 
orally, electronically, or even indirectly. Additionally, people interact accord¬ 
ing to different types of relationships, e.g. friendship and co-working, each 
one of these affecting the type of information transmitted. Consequently, 
a correct representation may require different types, or layers, of links [TB] . 
Neglecting such multi-layer structure, or in other words working with the 
projected network, may alter our perception of the topology and dynamics, 
leading to a wrong understanding of the properties of the system [HE]. 
In such cases, therefore, a single-layer network may be an oversimplihcation, 
and a multi-layer structure is required. 

A single-layer functional representation is generally created by collapsing 
the dynamics of the individual elements, such that the multiple dynamical 
aspects of each node are merged into a single time evolution. The same 
applies to the case of functional brain networks: the six-layer structure of the 
human cortex is neglected due to the limited spatial resolution of magnetic 
and electric sensors, and the analyzed time series correspond to the global 
activity of the top-most layers. Nevertheless, the non-linear nature of the 
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projection process can foster the appearance of constructive or destructive 
interferences. In other words, a link may appear in the projection even if no 
relationship is present in any layer; or links in two layers can interfere, and 
disappear from the projection. 

If one is to represent a real system by means of a functional network, 
and he / she expects it to have a multi-layer structure, it would then be 
important to understand the magnitude of the distortions created by the 
use of a single-layer representation. In other words, a fundamental question 
should be addressed: to what degree single-layer functional networks are 
representative of the dynamics occurring at different layers? 

In order to address this question, here I analyze the European air trans¬ 
portation network and create functional networks modeling the influence 
airports have on each another, in terms of flight delays and propagation. 
The availability of high-resolution real data allows the reconstruction of a 
complete multi-layer picture, in which each layer corresponds to a different 
airline; dynamics are afterwards collapsed, in order to simulate the creation 
of a single-layer representation. The resulting structures are compared, both 
topologically and dynamically, in order to assess whether they provide similar 
insight into the underlying system. 

2. Reconstructing the air transport functional network 

The air transportation system can be described both from a physical and 
a functional perspective. According to the former, the system supports a 
transport phenomenon, in which passengers and goods are transported be¬ 
tween pairs of airports. Therefore, the corresponding network representation 
encodes the presence of direct connections between nodes, in this case direct 
flights between airports. Numerous works have focused on this aspect, as for 
instance Refs. [iniiiTiii]. A latter perspective disregards the physical move¬ 
ment of items in the network, focusing instead on how the dynamics of nodes 
is influenced by other peers. Thus, a time evolving observable characterizes 
each element of the system, its dynamics being a function of the dynamics of 
its neighboring elements. When such functional relationships are uncovered, 
the result is a functional network. The air transport system can be seen as 
a collection of airports, whose observables are the delay experienced in each 
one of them. Clearly, delays are transported by aircraft, and thus the delay 
observed at one airport is a (potentially noisy) function of other airports de¬ 
lays. Such relationships are then susceptible of being mapped in a functional 
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complex network. 

The reconstruction of the physical network is a straightforward process, 
only requiring an analysis of flights’ schedules. On the other hand, the re¬ 
construction of delay functional networks starts with time series representing 
average landing delays across European airports. In this example, I have 
considered delay time series for the 50 busiest European airports and the 
20 largest airlines in number of flights during 10 consecutive months of year 
Appendix A for further details. These time series have been 
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pre-processed, in order to ensure their stationarity, for then assessing the 
synchronization level between pairs of them 


see 


Appendix B). 


As depicted in Fig. by starting from the raw time series (bottom left 
part), two possibilities can be followed. On one side, one can average all of 
the dynamics corresponding to a node, or in this case to an airport, thus 
creating a network representing the projection of the dynamics. Once one 
time series per airport is available, different metrics can then be calculated 
to construct a weighted fully-connected network. On the other side, one can 
directly create a multi-layer network by considering each airline as a layer, in 
which afterwards such structure can be projected into a single-layer graph, 
thus creating a projection of the topology. 

According to this idea, 23 different networks can be created: following 
the numeration of Fig. a projection of the dynamics (1); two projections of 
the topology, by respectively taking into account the strongest links in each 
layer (2) and the sum all statistically signihcant links (3), and the 20 airline 
networks (4 — 23). All networks, except for number 3, have been pruned in 
order to have a link density of 0.05. 


3. Topological properties 

Differences between the 23 reconstructed networks are hrst assessed from 
a topological point of view. Considering that this is a diffusion process involv¬ 
ing delays on top of the European air transportation network, the following 
topological metrics have been selected for their relevance: 

Out-degree centrality, i.e. the number of links departing from a given node. 
This centrality assesses the capacity of a node of inducing delays in 
neighbor airports. 

a-centrality, dehning the centrality of a node as the result of a random walk 
in the graph |18j . 


4 





I 


Projection of 
dynamics 



Airline delay time series 



LC 


GC 

TE 



Projection -> 2 3 

of topology 1 



Figure 1: Representation of the network reconstruction procedure. Starting from raw delay 
time series corresponding to different airports and airlines (bottom left), a multi-layer 
representation is created (bottom right). Furthermore, two projections are considered: 
a projection of the dynamics, where delays of airports are collapsed to create a single 
layer functional network (top left, and network number 1); and a topological projection 
of the multi-layer representation (networks number 2 and 3). In network 3, all links are 
considered, while in network 2 links are pruned to obtain the same link density as in 
network 1. LC, MI, GC and TF stand for the four synchronization metrics considered, as 
described in [Appendix B 


Efficiency, defined as the inverse of the harmonic average of the geodesic 
distance between all pairs of nodes, and measuring how efficiently the 
network can diffuse information [T9] . 

Clustering coefficient, a measure assessing the propensity of nodes to cluster 
together, and dehned as the number of closed triplets {i.e. the number 
of triangles in the network) over the total number of connected triplets. 

Degree-degree correlation, measuring the assortativity of the network, i.e. 
the fact that nodes tend to connect to their connectivity peers; it is 
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Table 1: Main topological properties of the reconstructed networks. From left to right, 
the six columns represent the following networks: the projection of the dynamics (network 
1 of Fig. [^, the two projections of the topology (2 and 3 of Fig. [^, and the first three 
layers (airlines). 



Proj. (1) 

Proj. (2) 

Proj. (3) 

Layer 1 

Layer 2 

Layer 3 

Link density 

0.05 

0.05 

0.54 

0.05 

0.05 

0.05 

Maximum out degree 

26 

16 

48 

13 

13 

15 

Degree-degree correlation 

- 0.058 

- 0.077 

- 0.029 

0.031 

- 0.088 

0.094 

Clustering Coefficient 

0.022 

0.179 

0.707 

0.064 

0.217 

0.095 

Efficiency 

0.064 

0.094 

0.765 

0.110 

0.064 

0.113 

Normalized IC 

0.875 

0.109 

0.179 

0.463 

0.298 

0.216 

Size giant component 

2 

17 

48 

16 

14 

18 


calculated through the Pearson correlation coefficient of the degrees at 
either ends of a link |20j . 

Entropy of the degree distribution, quantifying the network heterogeneity, and 
related with its robustness against random failures W- 

Normalized Information Content (IC), assessing the presence of regularities 
in the adjacency matrix, and thus the presence of meso-scale structures 
|22j . The lower the value, the more regular is the link arrangement, 
indicating the presence of meso-scale structures. 

Tableand Fig. [^report the main the results obtained in all considered 
networks. The two networks respectively obtained by projecting the dynam¬ 
ics and the topology of the system (number 1 and 2 of Fig. share a 90.5% 
of their links. Nevertheless, they are characterized by different topologies, 
as shown in Tab. [1} Specifically, the former has a centralized structure, as 
indicated by the higher maximum degree; the latter is less centralized, with 
a higher clustering coefficient. Due to the high number of links composing it, 
the network representing the unpruned projection of the topology (number 3 
of Fig. is the most different, sharing respectively only a 45.1% and 50.65% 
of links with the two peers. Such high connectiveness reflects in the topologi¬ 
cal metrics, yielding high clustering coefficient and efficiency. Finally, Tab. 
reports some topological characteristics of the hrst three layers, correspond¬ 
ing to the three most important airlines operating in Europe. While the first 
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and third one share most characteristics, the second one is characterized by 
an higher clnstering coefficient; this is dne to their different bnsiness models, 
being the first and the third major airlines, the second a low-cost. 

Figs. A and B focns on the distribntion of node centralities, thns 
indicating which nodes are most relevant from the point of view of delay 
propagation. In both panels, the main graph reports the centrality of the 
most important node in the network resnlting from the projection of the 
dynamics (network nnmber 1 in Fig. [^, as calcnlated by means of the ont- 
degree centrality (panel A) and a-centrality (panel B). The two horizontal 
lines and the vertical bars respectively represent the centrality of that node 
in the projection of the dynamics (by definition, eqnal to 1), in the projection 
of the topology (network 2 in Fig. [^, and in all the 20 layers of the mnlti- 
layer representation. Fnrthermore, the two small graphs on the right part 
depict the same information for the most central node in the first two layers. 
Figs. [^C, D and E report on three additional topological properties for the 
22 considered networks, i.e. efficiency, clnstering coefficient and entropy of 
the degree distribntion. 

In each case, a strong heterogeneity in resnlts is obtained. Specifically, 
the most important node in one network is seldom of major relevance in other 
networks. Also, topological features like efficiency and clustering coefficient 
take significantly different values depending on the type of projection con¬ 
sidered. This has important consequences from an engineering point of view. 
For instance, if one is to identify the most important node (airport) in terms 
of delay propagation, the use of different projections would yield different 
results, providing no guarantee that the most central node in the projected 
networks is indeed relevant when the multi-layer dynamics is considered. 

4. Dynamical properties 

From a dynamical perspective, it is necessary to ascertain whether such 
topological differences are enough to alter our understanding of the dynamics 
of the system. Specifically, in the case of air transport, this involves quanti¬ 
fying the amount of delay generated within the system, as well as its spatial 
and temporal evolution; on the other side, estimating the error resulting from 
the use of a projection of the multi-layer structure. 

In order to address these differences, a simple dynamical model has been 
constructed, mimicking the process of delay propagation in the real system. It 
comprises aircraft performing a random walk on top of the physical networks. 
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Figure 2: Characterization of the functional networks. (4) and {B) respectively repre¬ 
sent the out-degree centrality and a-centrality of the most central node in the dynamical 
projection (network 1 of Fig. in the network corresponding to the projected topology 
(network 2 of Fig. and in all layers of the multi-layer representation. The two small 
graphs on the right respectively represent the same results for the most central node for 
the first and second layer. (C), (D) and [E) respectively represent the efficiency, clustering 
coefficient and entropy of the degree distribution for the two projected functional networks 
and the layers of the multi-layer representation. 


with a stochastic delay added at each time step to simulate the appearance 
of uncorrelated disturbances. Aircraft delays are further modihed according 
to a delay multiplier | 23 ], which may take two different values: a when only 
a physical connection {i.e. a scheduled flight) is present between the start 
and destination nodes, and (3 when a functional link is also detected - see 
Appendix C| and Fig. |^A for further details. Thus a is associated to routes 


where the delay is expected to partly recover, while /3 indicates some kind 
of delay propagation. Two variants of the model are also considered. In the 
hrst one, all aircraft share the same network, i.e. the single-layer projection 
is used; in the second, the full multi-layer structure is implemented, allowing 
aircraft to respect the corresponding airline connection network. 

Once defined this two-fold model, it is relevant to study the error intro¬ 
duced by using the single layer representation, in terms of the amount of 
delay accumulated during the operation. In other words, we are interested 
in estimating the quantity Dml — Dsl (see Appendix C for the definition). 
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Figure 3: Construction and results of the dynamical model. {A) Nodes (airports) are 
connected both by physical links (black solid arrows) when aircraft are allowed to move 
between two nodes, and by functional links (green dashed arrow) when a functional re¬ 
lationship is detected between the corresponding delay time series. When both links are 
present, the delay multiplier has a value of /3, and a otherwise. [B) Difference between 
the delay obtained in the multi-layer and the single layer representations, as a function 
of a and /?. The two input parameters, a and /3, are respectively represented by the X 
axis and by different curves. To obtain the estimation error, one thus has to fix the X 
position (a), intercept the desired curve (/3), and read the result in the V axis. (C, D) 
Probability distribution for links in the single and multi-layer networks; light and dark 
green respectively represent physical and functional links. 


and in describing how it evolves as a fnnction of the two parameters a and 
(3. The resnlting estimation error is presented in Fig. B as a fnnction 
of a and (3. For large valnes of l3, the delay obtained in the mnlti-layer 
representation increases to twice the value obtained in the single-layer case. 
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Figure 4: Asymptotic dynamics of the dynamical model. Panels A and B respectively 
represent In 7 as a function of a and /3, for the single layer and multi-layer models. For 
the sake of clarity, values of 7 smaller than 1 have been set to 1 , such that ln 7 > 0 . 


Such error is due to the different topology of both networks: specifically, a 
significant number of functional links [i.e. associated to the /3 multiplier) in 
the multi-layer representation have a high centrality, thus indicating that key 
airline connections are also the ones responsible for the delay propagation; 
on the other hand, such connections may loose importance in the projected 
networks, thus reducing the amount of delay generated (see Fig. |^C and D). 

In order to better understand the differences between the single and multi¬ 
layer scenarios, the average increment of delays (denoted by 7 ) has been 
Appendix C Fig. [^reports the value of In 7 as a function of 


see 


calculated 

a and {3 - the two delay multipliers, as in Fig. The single and multi-layer 
scenarios (respectively left and right panels) present two distinct behaviors. 
Specifically, the change from a steady state (7 ~ 0) to a regime of increasing 
delays (7 > 0 ) depends both on a and f3 when the dynamical model is 
executed on a single-layer projection. On the contrary, the behavior of the 
multi-layer model strongly depends on the value of a, (3 having a neglectable 
impact. 


5. Discussion 

While functional networks have represented a great leap forward in the 
understanding of complex systems, the most notable of them being the hu- 
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man brain, little attention has been devoted to the effects of collapsing multi¬ 
layer dynamics into a single layer functional graph. When the dynamics of 
different layers have similar characteristics, the use of a projection is ex¬ 
pected to have negligible effects on the results; nevertheless, when layers are 
heterogeneous in nature, such simplification may be misleading. 

In this paper, 1 tackle this issue by analyzing the problem of delay prop¬ 
agation in the European air transport network. Different functional network 
representations are compared, i.e. the complete multi-layer representation, 
and the corresponding projections of the topology and dynamics. When an¬ 
alyzing their topological characteristics, important differences are observed. 
Notably, the most central nodes in the projections do not correspond to the 
nodes of high centrality in each layer; therefore, the former analysis give 
biased estimations, which cannot reliably be used to detect the most criti¬ 
cal elements in the system. Furthermore, when a simple dynamical model is 
executed, the magnitude of the error yielded by considering a single layer pro¬ 
jection is as big as the results themselves, thus indicating that any estimate 
obtained with this simplification is meaningless. 

These results present important implications for the modeling and fore¬ 
casting of the air transport network. First, they imply that any simulation 
performed to understand the dynamics of the system may yield misleading 
results when the multi-layer structure created by airlines is neglected. In 
spite of this, most of the recent research works in this fields fail to include 
this essential ingredient, both in the analysis of delay propagation [231 ED 
and of the network robustness to disruption and attacks [25l|26l|27|. Second, 
it has to be noticed that the air transport system is created by the inter¬ 
actions between a large number of agents, which may create different layers 
along different dimensions. For instance, multiple flights do not just share 
the airline, but they may also be connected by the crew operating them. 
Disregarding these different layer dimensions, e.g. crews, aircraft types or 
flight type (cargo or passengers), may further bias our understanding of the 
system. While it may not be possible to include all these elements in one 
analysis, for instance due to the lack of publicly available data, it is safe to 
conclude that “one layer does not fit all”. 

Beyond the air transport network, the results here reported can be ap¬ 
plied to any other multi-layer complex system analyzed by means of func¬ 
tional networks: the most paradigmatic example is then the study of the 
brain. While the air transport network and the human brain are different in 
nature, some results obtained for the former can be translated into the anal- 
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ysis of the latter. Specifically, both systems are characterized by a transport 
layer in which items are transported between nodes. The structure created 
by direct flights is conceptually similar to the connectome, i. e. the structure 
created by hbers connecting different regions of the brain, which allows in¬ 
formation (spike trains) to move between them [2H]- Both systems can also 
be analyzed from an information processing point of view: the air trans¬ 
port processes information about delays, as transported by flights, while the 
brain information about the body and the environment, codihed by neural 
spikes. If the most important airports, in term of delay propagation, cannot 
reliably be detected with a projected functional network, the identihcation 
of functional hubs in the brain dynamics may be confused by the fact that 
the multi-layer structure of the cortex is neglected. Therefore, global hubs 
may not correspond to the most important nodes in each layer: the single 
layer analysis may then be misinforming about the real structure created by 
information flows. In order to conhrm and complete our knowledge about 
brain dynamics, it is then desirable to develop techniques able to accurately 
represent its multi-layer structure and dynamics. 
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Appendix A. Air traffic data 

Information about aircraft trajectories has been extracted from the Flight 
Trajectory (ALL-FT-I-) data set as provided by the EUROCONTROL PRISME 
group. It includes information about planned and executed trajectories for 
all flights crossing the European airspace, with an average resolution of 2 
minutes. The data set covers the period from 1'^* March to the 31®* De¬ 
cember 2011, including a total of 10.3 million flights. From this, a subset 
has been considered, composed of the 50 busiest airports and the 20 largest 
airlines (in number of flights) operating at these airports. 

Appendix B. Time series extraction 

For each combination of airport / airline, a time series has been extracted, 
representing the average hourly delay of arriving flights. Due to some miss- 
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ing days, each time series comprises 7440 values. The resulting time series 
presented several trends, as strong delays are expected mainly at peak hours, 
during week days, and during summer, i.e. those periods in which the traffic 
is higher. A detrend process has then been performed, by subtracting the 
average delay observed in the same day, in the two previous and posterior 
weeks, at the same hour, he.: 

d{t) = d{t) - ^ E d(t + 168i), (B.l) 

d{t) being the original time series at time t, and d{t) the hnal time series. 
According to this dehnition, d{t) thus represents the difference between the 
observed and the expected (historical) delay. 

Functional networks have been reconstructed by considering all pairs of 
airports, and by assessing the presence of a relationship among their time 
series by means of four metrics: 

Pearson’s linear correlation (LC): a well-known measure of linear correla¬ 
tion, is dehned as the covariance of the two variables divided by the 
product of their standard deviations [2U] . 

Mutual Information (MI): measures the information shared by two random 
variables, X and Y, by analyzing their individual and joint probability 
distribution functions [50] : 


= (B.2) 

\p(Op(.y)J 

Granger Causality (GC): measures the presence of a causal relation, i.e. if 
one time series is useful in forecasting a second one [I2] . Given two time 
series X and T, X causes Y ii X values provide statistically signihcant 
information about future values of Y. This is usually calculated by 
means of univariate autoregressions, and tested through t-tests and 
F-tests. 

Transfer Entropy (TE): given two processes (or time series) X and Y, TE 
is the amount of uncertainty reduced in future values of Y by knowing 
the past values of X and Y: 
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H being the Shannon’s entropy of the time series. While similar in 
nature, the advantage of TE with respect to GC is that the former can 
be used in the analysis of non-linear signals inn. 

These four synchronization metrics have been chosen to cover different 
characteristics of the time series: they include two correlations (LC and MI) 
and two causalities (GC and TE); furthermore, two of them are linear (LC 
and GC), two non-linear (MI and TE). Results obtained with these four 
metrics are qualitatively similar, thus only the ones corresponding to TE are 
here reported. 

Finally, a physical network has been reconstructed, where nodes represent 
airports, pairs of them being connected by unweighted links when at least one 
direct flight has been operated between them. Two variations of this network 
are considered for constructing the dynamical model (see next Section): the 
complete multi-layer structure, each layer representing the network of flights 
of a single airline, and a single layer projected network. 

An interactive representation of the obtained networks, both physical and 
functional, can be found at |32j . 

Appendix C. Dynamical model 

In the dynamical model representing delay propagation through the air 
transportation system, a set of n = 10.000 aircraft are randomly placed in 
the 50 nodes of the network; in the case of the multi-layer representation, 
aircraft are also randomly assigned to one of the 20 layers available. Initially, 
all aircraft are on time, i.e. their delays di{t = 0) = 0 with i G [1,^.]. At 
each time step t > 0, aircraft are randomly moved according to the physical 
adjacency matrix, corresponding to the flight connections detected in the 
system, respecting the assigned layer in the case of the multi-layer structure; 
furthermore, their delays evolve according to: 


di(t) = adiit — 1) -|- e, if a = 0, 
diit) = jddiit — 1) -|- e, if a = 1, 


diit) = 0,if diit) < 0, 
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(G.l) 

(G.2) 

(G.3) 


e being a random number drawn from a normal distribution 7V(0,1), and 
a representing the presence of a functional link between the start and desti¬ 
nation nodes (see Fig. [^A). The last condition simulates the scheduling of 
flights: when an aircraft reaches its destination airport before the scheduled 
time, it waits until the expected time of departure. 

The two parameters a and (3 dehne how delays are expected to evolve 
after every flight, and are commonly known in the aeronautical community 
as delay multipliers [23]. A value larger (smaller) than 1.0 indicates that each 
flight is, on average, arriving with a delay that is larger (respectively, smaller) 
than the one recorded at take-off, and is thus multiplying (dampening) the 
total delay. In order to simulate this process in the dynamical model, each 
possible route is associated to either a or /?, depending on the existence of the 
corresponding functional link between the departure and destination airports 
(see Eq. C.2 and C.3). When a = 0, i.e. no causality link has been detected 


between two airports, the delay propagation is controlled by a; the lack of 
a causality relation suggests that delays are (at least partly) recovered, and 
thus a is expected to be smaller than 1.0. On the other hand, if a = 1, a 
delay propagation process is probably present between the two airports: (3 
is thus expected to be greater than 1.0, and also that f3 > a. It should be 
noticed that the single and multi-layer networks may have a different number 
of physical and functional links. In order to have comparable results, in the 
latter case the model is executed with two parameters a* and 13* dehned as: 


/3* 

p and / respectively being the proportion of physical and functional links 
in the multi-layer [ML) and single layer networks {SL). 

The hnal delay generated by the system is measure as: 

D = ^d.(t = 100). (C.6) 

i 

Results reported in Fig. B correspond to the average total delay D 
obtained in 100 realizations of the model. 

The importance of each link in the networks (Fig. [^C and D) is calculated 
as the a-centrality of the corresponding dual graph, where link weights are 






= /9 


fML/fSL 


(C.4) 

(C.5) 
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normalized such that the out-strength of each node is equal to one - thus 
simulating a random walk on the network. 

Finally, the average increment of delays has been calculated as: 

(C.7) 

i i 

di(t) being the delay observed at time t for airport i, and (■) the average 
through time. Values of 7 close to 0 imply that the system is able to dissipate 
the delays generated by the random part of Eqs. 4-6; on the contrary, for 
7 > 0 , a steady state is not reached and delays constantly increase. 
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